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Abstract: We study the velocity-dependent cusp anomalous dimension in supersymmetric 
Yang-Mills theory In a paper by Correa, Maldacena, Sever, and one of the present authors, 
a scaling limit was identified in which the ladder diagrams are dominant and are mapped 
onto a Schrodinger problem. We show how to solve the latter in perturbation theory and 
provide an algorithm to compute the solution at any loop order. The answer is written in 
terms of harmonic poly logarithms. Moreover, we give evidence for two curious properties 
of the result. Firstly, we observe that the result can be written using a subset of harmonic 
polylogarithms only, at least up to five loops. Secondly, we show that in a light-like limit, 
only single zeta values appear in the asymptotic expansion, at least up to six loops. We then 
extend the analysis of the scaling limit to systematically include subleading terms. This 
leads to a Schrodinger-type equation, but with an inhomogeneous term. We show how its 
solution can be computed in perturbation theory, in a way similar to the leading order case. 
Finally, we analyze the strong coupling limit of these subleading contributions and compare 
them to the string theory answer. We find agreement between the two calculations. 
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1. Introduction 

The cusp anomalous dimension r cusp (</>) was originally introduced in [1] as the ultraviolet 
(UV) divergence of a Wilson loop with a cusp with Euclidean angle 4>. It describes a wide 
range of interesting physical situations. It was computed in QCD to the two-loop order in 
ref. [2] and rederived and simplified in ref. [3]. 

In super symmetric Yang- Mills theories such as N = 4 super Yang-Mills, one can define 
a Wilson loop operator that couples to scalars in addition to the gauge field [4,5]. It is 
natural to consider a loop where the coupling to the scalars is different on the two segments 
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of the cusp (but constant along each segment). The jump in the internal coupling to the 
scalars is characterized by an angle 9. The perturbative calculation of this super symmetric 
cusp anomalous dimension T cusp (4>,9) is similar to the QCD case. It has been performed 
to two loops in refs. [6,7]. At strong coupling, it is known to second order in the strong 
coupling expansion [7]. 

Recently, there has been a lot of progress in understanding r cusp (0, 9), in various 
domains. 

In a small angle angle limit, an exact result was found in [8], based on localization 
techniques. The exact formula is in perfect agreement with perturbative results and the 
result at strong coupling. The same exact formula has also been obtained in [9] . 

In ref. [10], a relation of the cusp anomalous dimension to the Regge limit of mas- 
sive scattering amplitudes was exploited to compute its three-loop value. The relation to 
scattering amplitudes [11], which is valid in the planar limit, implies in particular that 
the integrand needed to compute the cusp anomalous dimension can be deduced from the 
(in principle known) integrand for planar four-particle scattering amplitudes [12-14], when 
appropriately extended to the massive case [11,15]. 

Very recently, Thermodynamic Bethe Ansatz (TBA) equations have been derived for 
the cusp anomalous dimensions [16,17], and passed highly non-trivial consistency checks 
at the three- loop level [16]. 

In [10] a new scaling limit involving the complexified angle 9 was introduced, 

i9 > 1, A < 1 , with A = Xe ie /4 finite. (1.1) 

Here A = g 2 N is the 't Hooft coupling. In this limit, the coupling of the loop to the 
scalars becomes dominant, and the leading order (LO) contribution is given by simple 
ladder diagrams, where the rungs of the ladder are scalar exchanges. It is important to 
realize that this is a gauge-invariant statement. The ladder diagrams can be described 
conveniently using Bethe-Salpeter equations. The latter are very convenient, since they 
provide a simple description. They can be solved exactly in the small angle limit, and it is 
easy to extract their strong coupling behaviour, finding agreement with the corresponding 
string theory calculation. 

In this paper, we continue the analysis of the scaling limit of [10] and initiate a system- 
atic study of the subleading contributions. A first question that one faces when computing 
rcusp(0, 6) in perturbation theory is what functions the result can be expressed in. It is 
easy to see that the 6 dependence is very simple, and to describe the <p dependence the 
variable x = e 1 ^ is useful. Experience shows that in that variable one obtains certain 
polylogarithms, multiplied by rational prefactors. In general, it is not known what class of 
polylogarithms, or more generally what class of iterated integrals, is sufficient to describe 
a given problem. 

Similar questions are of great current interest in the understanding of the structure 
of scattering amplitudes, a problem that is closely related. To phrase the question in 
that language, given a loop integral depending on n space-time points, what is the set of 
functions describing it? On the one hand, one could argue that with increasing loop order, 
integrals with some number of external points "know about" lower-loop integrals with 
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more external points that they contain as subdiagrams (which may e.g. contain elliptic 
integrals), making them very complicated, and perhaps requiring a larger functional basis 
at higher loops. On the other hand, one might argue that the set of functions should 
ultimately be determined by the external kinematics of the problem. An argument in favor 
of this point of view is that integrals are determined to a great part by their singularities, 
and the location of the latter is intimately tied to the external data. These questions are 
also of enormous practical importance, as they sometimes allow to make an ansatz for a 
given problem within a restricted class of functions, see [18-20] for recent examples. 

In the present case of a single scale problem, it was observed in [10] that all functions 
occurring to the three-loop order could be expressed in terms of harmonic polylogarithms, 
i.e. in terms of iterated integrals with integration kernels 1/x, l/(l + x), 1/(1 — x). The fact 
that this was possible not only for the final answer, but also for individual loop integrals, 
and in fact also for all integrals of this type found in the literature, seems to suggest that 
this is a more general feature. Can this be proven rigorously? In this paper, we make a 
first step into this direction. We show that this property holds for the LO term of r cusp in 
the scaling limit (1.1), and for one of the two contributions at NLO, at any loop order. 

We also present an algorithm that determines r cusp at LO in the scaling limit at any 
loop order in terms of harmonic polylogarithms. As an application, we verified the result 
of [10] at three loops, and evaluated the four- and five-loop results, which are new. 

These results suggest two further properties. First, we find that at least up to five loops, 
one can express the result in terms of harmonic polylogarithms (HPLs) [21] of argument x 2 
and indices 0, 1 only. Second, in the x — > limit we find that, at least up to six loops, single 
zeta values and products thereof are sufficient to describe the coefficients of the asymptotic 
expansion. 

We then discuss NLO terms in the scaling limit. We show that there are two classes 
of diagrams that satisfy a slightly modified Bethe-Salpeter equation. For one of the two 
classes of integrals, we show how to construct the solution in terms of HPLs at any loop 
order. For the second class of integrals, we compute the non-trivial integration kernel, 
which allows to express the result in terms of iterated integrals having the correct degree. 
We leave the question of whether the latter can be expressed in terms of HPLs to future 
work. 

We also discuss the strong coupling limit of the Bethe-Salpeter equations, and compute 
the scaling limit of the corresponding string theory result. Under certain assumptions, we 
find perfect agreement between the two calculations. 

This paper is organized as follows. We begin by reviewing the definition of the cusp 
anomalous dimension and the scaling limit in section 2. Then, in section 3, we present 
the perturbative solution at leading order in the scaling limit to any loop order. We prove 
that the result can be written in terms of HPLs, and make further observations about their 
structure. In section 4, we discuss the NLO Bethe-Salpeter equations. In section 5, we 
take the strong coupling limit of the equations, and compare them to the corresponding 
string theory calculation. There are several appendices containing technical details. 
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(a) (b) (c) 



Figure 1: Classes of loop integrals contribution to LO (diagram (a)) and NLO (diagram (b) and 
(c)) in the scaling limit ([O]). Each class can have an arbitrary number of rungs. The arrows in (c) 
denote a numerator factor (p + q) 2 , where and q^ are the momenta along the arrows. 

2. General structure of the Bethe-Salpeter equations at LO and NLO 

Here we discuss the general structure of the Bethe-Salpeter equations at leading order 
(LO) and next-to-leading order (NLO) in the scaling limit. The LO equations were already 
discussed in [10]. They are a natural generalization of the equations for the quark-antiquark 
potential [22]. Here we briefly review the main points. 

We recall the definition of the locally super symmetric Wilson loop operator in N = 4 
super Yang-Mills, 

W ~ ^p e i § A-dx+§\dx\n-^ ^ £2.1) 

where n is a point on S 5 . The contour we consider consists of two (infinite) segments 
forming a cusp of Euclidean angle <f>. We take the coupling to the scalars to be constant 
along each segment, but with a jump of angle 8 at the cusp, i.e. cos 9 = n ■ n' , where n 
and n! are the directions of the two segments. Such a cusped Wilson loop in general has a 
logarithmic divergence that takes the form 

(P^) ~e" rcusp(</,,9)log ^, (2.2) 



where Ajr/uv are infrared and ultraviolet cutoffs, respectively. This defines the cusp 
anomalous dimension r cusp (0, 9). 1 

In the scaling limit (1.1), the scalar coupling of the loop becomes dominant. At leading 



order (LO) in the limit, the segments of the Wilson loop couple to conjugate scalars, and 
we need to consider scalar exchange diagrams only. At next-to- leading order (NLO), we 
have mostly scalar exchanges, plus one-loop interaction diagrams. 

An analysis of the integrals contributing to the cusp anomalous dimension allows one 
to see that the effective diagrams shown in Fig. |l| are needed at LO and NLO in the 
scaling limit. Since only one-loop internal graphs are allowed at NLO order, one can 



x Of course, F CUBp is also a function of the 't Hooft coupling g 2 N, and the number of colours N. 
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Figure 2: Bethe-Salpeter equation at LO and NLO. The arrows denote a numerator factor (p + q) 2 , 
with being the momenta flowing along the arrows (in momentum space). 

deduce the all-loop structure of these corrections already from the known three-loop ex- 
pression. The fact that one has effective diagrams that arise after cancellations between 
various gauge-dependent Feynman diagrams 2 is intimately related to the similar diagrams 
appearing in scattering amplitudes. We illustrate this relation at the level of the loop 
integrals/integrands in Appendix C. 

It is easy to see that the integrals of Fig. [I] are described by a Bethe-Salpeter equation. 
The latter is shown (schematically) in Fig. ||. This equation sums the diagrams to all orders 
in the coupling. At LO in the scaling limit, only the first line contributes, as the second 
lines gives contributions of order a = A/A and higher. At NLO, we keep the terms in the 
second line and compute the answer linear in a. Note that there are also higher-order 
terms in a contained in this equation that will only become relevant once we include all 
NNLO and higher terms. 

We can see that there are two new features w.r.t. LO. First, the first term of the 
second line of Fig. ^ is the starting point for the new infinite class of diagrams shown in 
Fig. 0(b). These terms are absent in the quark-antiquark potential [23]. Second, there is a 
new interaction term that is a higher-loop generalization of the simple scalar exchange at 
LO. 

Let us illustrate the usefulness of the Bethe-Salpether equation by reviewing the LO 
case. We denote the sum of the ladder diagrams by F(s,t), where — sp^ 1 and tq^ are 
positions on the cusp formed by the momenta p^ and q^. Let us normalize p 2 = q 2 = 1 for 

2 In ref. [23], this one-loop calculation was explicitly performed (for the quark-antiquark potential, cor- 
responding to <f> — > 7r), in agreement with the result here. Integral class (b) discussed here follows from a 
boundary term at the cusp that is absent for the quark-antiquark potential. 
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convenience. Note that F also depends on the angle <f> defined by cos</> = p ■ q. Then F 
satisfies the Bethe-Salpeter equation 

F(S,T) = 1+ [ ds [ dtF(s,t)P(s,t) , (2.3) 
Jo Jo 



where 



P(s '* )_ 47r 2 S 2 + t 2 + 2 S tcos</> (2 ' 4) 

is the propagator corresponding to a scalar exchange. Changing variables according to 
s = e a , t = e T , this becomes 

/cr f'T 
da 1 / dnFiaur^P^n), (2.5) 
-oo J — oo 

where 

^^A-i ^— -■ (2.6) 

07r z cosn(r — a) + cos 

Differentiating eq. (^Jj), we obtain, 

d T d a F(a, t) = F(a, r)P(a, r) . (2.7) 

Let us change variables yi = r — a and yz = (r + <r)/2. We can extract r cusp from the 
large y2 behaviour of F, due to the equivalence of IR and UV divergences, see eq. (f2^). 
For large y2, we can make an ansatz 

F = Y,e- n " y ^n(yi)- (2.8) 

n 

We are interested in the leading term, corresponding to the lowest eigen-energy Qq. Using 
the ansatz (2J3), one finds [10] 

*(yi,^) = 0. (2.9) 



2 A 1 O 2 (0) 



8ir 2 (cosh yi + cos i 



This is a one-dimensional Schrodinger problem. The ground state energy is related to 
the cusp anomalous dimension in the scaling limit through r cuS p — 

In summary, the Bethe-Salpeter equation has allowed us to conveniently sum an infinite 
class of diagrams. As a result, extracting the remaining overall logarithmic divergence could 
be done in a simple way, and the remaining calculation does not require any regulator. 
Moreover, the structure of the equation allowed us to rewrite the problem in terms of a 
linear differential equation. 

We will now solve this equation in perturbation theory. In section |], we will discuss 
the effects of the two new features that appear at NLO. 
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3. Solution to the scaling limit at leading order 



3.1 Setup 

To obtain the perturbative solution of (|2.9|) , we follow [10] and perform the change of 
variables 



(3.1) 



The exponential factor gives the correct solution as y\ — ¥ oo, and we can normalize 77(7/1 
00) = 1. We can determine VLq from r] thanks to the boundary condition 



dyMv 



U 1 3/1=0 



0. 



(3.2) 



which follows from the y\ — > —y\ symmetry of the problem. Defining a new variable 
w = e~ Vl , and x = e"^, the boundary condition ( |3.2| ) becomes 



£l (x) = -2wd w log 7](w, x)\ w =i , 
and the Schrodinger equation ( [2.9D reads 

1 1 



d w wd w rj = -£l (x)d w r] + k 



Xx 



V: 



K 



47T 2 (1 - X 2 ) 



(3.3) 



(3.4) 



W + X -1 W + X 

The wavefunction r/ can be obtained by integrating the Schrodinger equation iteratively in 
l -r n -r . . ., aii>a </ — x t rv</ • • t Let us now analyze in 



the coupling, Qq = kQ,^ + k 2 Q,^ + . . ., and r\ = 1 + krf^> + 



detail the perturbative solution for rj and 0. 
3.2 Iterative solution 

It is convenient to introduce an abbreviation for the nested integrals that one encounters 
in this problem. In analogy to two-dimensional harmonic polylogarithms (2dHPLs), we 
are going to use the self-explanatory notation 

Hv(w, x) 



dw' 



1 



1 



and 



H v r( W ,x) 



dw' 



w' + l/(wx) w' + x/w 

HAw'w, x) 



w' + l/(wx) w' + x/w 
f 1 dw' 



(3.5) 

(3.6) 
(3.7) 



In the following we will sometimes drop the arguments (w, x) for brevity. So in general we 
will have H^, where the weight vector b has entries V and 0, with not appearing in the 
last entry. 

It is straightforward to write the perturbative answer for r\ in terms of these integrals. 
We find 



=H 0jV 

V {2) =H ,v,o,v - H , , v nP , 



(3.8) 
(3.9) 
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and so on. Using eq. (^) we find 

n { 1] =-2H v , (3.10) 

n { Q ] = - 2H v H Q ,v - 2H VfitV , (3.11) 
Oq 3 - 1 = — 2HvHq V — 4HyH 0: oy — 2H yHvfiy 

— 2iJyHo,V,0,V — 4HyHv,0,0,V — 2Hv,o,v,o,v , (3-12) 

etc. These last relations are understood at w = 1. 

In principle, eqs. ( |3.10|) , (|3.1l|) , ( |3.12|) , and their higher-order analogues, together with 
r = —^0) provide formulas for T. However, this representation is clearly not an optimal 
one. In the following, we will simplify it by converting it to a more appropriate and simpler 
class of iterated integrals. This will also allow us to make further observations regarding 
the structure of the result. 

3.3 Structure of the perturbative result 

Here, we first show certain properties of r/ and Jlo, and then outline an algorithm for 
expressing Q$ in terms of harmonic polylogarithms. 

As we show presently, the total differential of rj at any loop order is of the form 

dr/ L ) =f l d\ogx + f 2 dlog(l + x) + h dlog(l - x) 

+ hd\o%{w + x) + fodlog{w + l/x) , (3.13) 

with the fi being functions of the same type as r/ L ) , but of degree (i.e. number of iterated 
integrals) lowered by one. From equations ( 3.13j ) and ( |3.3[ ) it immediately follows that 



dn^ = gi d\ogx + g 2 d\og(l + x) + g 3 dlog(l - x) , (3.14) 

with gi being functions of degree lowered by one, and satisfying the same property. This, 
implies that at any loop order L, fj( L ) can be expressed in terms of harmonic polylogarithms 
(HPLs) of degree (2L-1). 

The latter are defined iteratively by 

H ai ,a 2) ...,an(x) = I f ai (t)H a2> ... >an (t)dt, (3.15) 
J 

where the integration kernels are 

/ 1 (x) = -J_, fo(x) = -, f_ 1 ( x ) = J—. (3.16) 

1 — X X 1 + X 

The degree-one functions needed to start the recursion are defined as 

H 1 {x) = -\og(\-x), H (x) = log(x) , ff_i(x)=log(l + x). (3.17) 

The subscript of H is called the weight vector. A common abbreviation is to replace 
occurrences of m zeros to the left of ±1 by ±(m+l). For example, ffo,o,i,o,-i( a; ) = ^3,-2(^)- 
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Note that a corollary of equations ( 3.13 ) and ( 3.14 ) is that the symbol [24,25] of rj is 
constructed from a five-letter alphabet consisting of x, 1 ± x, w + x, w + 1/x. Similarly, the 
symbol of f^o is constructed from the three-letter alphabet x, 1 ± x. Of course, knowing 
the full differential provides us with much more information than just the symbol. 

In order to prove the above statements, let us point out a relation of the H^(w, x) to 
a known, albeit more general class of functions, the Goncharov polylogarithms [26], 



G(ai, . . . a n ; z) = 

Jo 



dt 



t — a\ 



G(<22, . . . , a n ; t) , 



with 



G(ai;z) 



dt 



io t — a\ 

In our case, the are taken from {0, —x, —1/x} and z = w. For example, we have 

Hy(w, x) =G(— 1/x; w) — G(—x; w) , 
Hoy(w, x) =G(0, —1/x; w) — G(0, — x; w) , 
Hvfi,v(w, x) =G{—l/x, 0, — 1/x; w) — G{— 1/x, 0, — x; w) 
+ G{— x, 0, — x; w) — G(—x, Q,—l/x;w), 

and so on. The total differential of a general Goncharov polylogarithm is 
dG(ai, . . .a n ;z) = G(ai, a 2) - - - a n ; z) cZlog "' 



(3.18) 



(3.19) 



(3.20) 
(3.21) 

(3.22) 



a\ — 02 

+ G(ai, a 2 , a 3 , . . . , a n ; z) (flog — — ^ + . . 



+ G(oi, . . . ,a n _i,a n ; z) dlog 



a 2 - a 3 



(3.23) 



where a means that this element is omitted. 

Given the possible values of the dj in our case, it is straightforward to verify eq. ( 3.13j ). 



3.4 Rewriting the expressions for f^o in terms of HPLs 

We have proven that f?o can be written in terms of HPLs. Let us now explain how to find 
explicit results in terms of HPLs. We will begin by a simple example, and then outline an 
algorithm for doing so in general. 

We observed that eq. ( p. 23 ), when applied to any function Hr(w = 1, x) gives a result of 
the form ( 3.14| ). Iterating this procedure for the lower degree functions gi in that equation, 
together with the fact that at any order we have a boundary condition at w = 1, gives us 
the complete information for that function, in a form that makes contact with the definition 
of HPLs, see eq. (3.15). 

As an example, let us write Hoy(l,x) in terms of HPLs. According to eq. ( 3.21| ), we 
need to rewrite G(0, — x; 1) and G(0, —1/x; 1) in terms of HPLs. Specializing ( 3.23|) to the 
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present case we have 

dG(Q,-x;l) = -G(-x;l)dlogx 

= — log((l + x)/x) (flog x 

= - [fT_i(x) - H (x)) dlogx. (3.24) 
The integration can be done using the definition ( 3.15| ), 

G(0, -x- 1) = - H ,-!(x) + H , (x) + C . (3.25) 
It is convenient to relate G to the value at x = 1, 

G(0, -x; 1) = - H Qi - X {x) + H 0>0 (x) + ^2 + G(0, -1; 1) , (3.26) 

where we used that i?o,-i(l) = V^C2- 
Similarly, we find 

G(0, -1/x; 1) =H ^(x) - ^ 2 + G(0, -1; 1) . (3.27) 
Combining eqs. ( |3.26| ) and ( p. 27 ) we find 

H oy (l,x) = 2H ,-i(ar) - #o,o(*) - C2 • (3.28) 

By construction, Hoy(l, 1) = 0. 

In summary, from this example it becomes clear how to rewrite any of the functions 
occurring in our problem in terms of HPLs, using the following steps: 

1. Express the functions in terms of Goncharov polylogarithms 

2. Use eq. ( [3.23 ) in order to compute their (total) differential; since all other variables 
are constants, this gives the derivative in x. 



3. By iteration, that differential is of the form (3.14), with the gi appearing there being 
HPLs. It can therefore be integrated in terms of HPLs, using (3.15). 

4. Integrate the equation with the boundary term at x = 1. 

5. Add up all terms; the boundary Goncharov polylogarithms at x = 1 do not necessarily 
drop out, but they are simple since they correspond to harmonic polylogarithms 
evaluated at x = 1. 

Using this algorithm, we find e.g. the following expressions that are required to the 
three-loop order, 

H v (l,x)=H (x), (3.29) 
H VAV (l,x) = - 4F_ 3 (x) - (2H (x) + 2#_ 2i0 (x) - 4H 2)0 {x) - H 0A0 (x) - 2( 3 , (3.30) 

H ,o, ,v(h x )=~^ + 2#-4(a0 - C2#o,o(x) " H ,o AO (x) , (3.31) 
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H ,v,o,v(l,x) - 14#_ 4 (x) - tt 2 H^ 2 (x) + \ir 2 H 2 {x) + 8#_ 3 ,-i(x) - 4#_ 3)0 (x) 

360 o 

+ 12#_ 2 _ 2 (x) + ^vr 2 ^ ,o(x) - 8# 2 ,_ 2 (x) - 6ff_ 2 , , (x) 
o 

+ 4F 2 , ,o(x) + #0,0,0,0^) + 2H Q (x)( 3 , (3.32) 

2tt 2 197T 4 47T 2 

# WA y(l,z) =40#_ 5 (x) - — #- 3 (x) + -5^H (x) + —H 3 (x) + 24#_ 4 ,_ 1 (x) 

o ooO o 

- 38#- 4 ,o(z) - 16#_ 3 ,_ 2 (x) - 24#_ 2i _ 3 (x) - ^#-2,0(2) + 32# 2i _ 3 (x) 
4tt 2 

+ — H 2 , (x) + 32# 3> _ 2 (x) + 52# 4 , (x) + 8#_ 3 _i )0 (x) - 4#_ 3 ,o,o(x) 

7T 2 

- 16#_ 3 , lj0 (x) + 12#_ 2 ,_ 2i0 (x) - 24F_ 2j2i0 (x) + y#o,o,o(x) 

- 16# 2 ,_ 2i0 (x) + 32# 2i2i0 (x) - 16H 3 _ li0 (x) + 8# 3A0 (x) + 32H 3jl>0 (x) 

7r 2 C 3 

- 6H_2,o,o,oW + 8H 2: oflfl(x) + #0,0,0,0,0(2;) + — 12C 3 i?- 2 (x) 

+ l6( 3 H 2 (x) + 2C 3 F ,o(x) + 6C5 ■ (3.33) 

Plugging these formulas into eq. ( |3.12D , we find perfect agreement with the three-loop 
result of ref. [10] . 

In the next section, we show explicit new results that we obtained using this algorithm. 



Using the method described in the previous section, we explicitly determined Q^\x) 



3.5 Explicit new results, and further surprises 

Jsin 

fig (x) in terms of HPLs. We will show these formulas below. 

When analyzing the resulting formulas, in fact we found a further simplification, that 
was already noticed in [10] up to the three-loop level. Although results for individual 
integrals contain in general HPLs with all possible indices 0, ±1, we observe that, at least 
up to five loops, it is possible to write the final result in terms of HPLs having indices 0, 1 
only, provided that we use x 2 as argument instead of x. That property is manifest in the 
following formulas. For completeness, we give all functions through to five-loop accuracy. 
Up to three loops, one finds 

^ 1} (x)= -# 0) (3.34) 
n { 2) (x) = 4 Cs + 2 C2 #o + 2 F 2 ,o + #0,0,0 , (3.35) 
n { 3) (x)= - 8 C2 Cs - 12 Cs - 12 (4 Ho - 16 Cs #2 - 8 C2 #3 - 4 Cs #0,0 - 8 C2 #2,0 

— 8 H^o — 8 C2 #0,0,0 — 8 #2,2,0 — 4 # 3 ,o,o — 8 #3,1,0 — 4 #2,0,0,0 — 6 #0,0,0,0,0 • 

(3.36) 

Our result at four loops reads 

S^ 4) (x) = 48 Cs + 24 C2 C 5 + 36 ( 7 + 8 C 3 2 #0 + 51 Ce #0 + 48 C2 Cs #2 + 72 Cs #2 
+ 96 C 4 #3 + 88 Cs #4 + 80 C2 #5 + 32 C2 Cs #0,0 + 20 C 5 #0,0 + 72 C 4 #2,0 
+ 96 Cs #2,2 + 48 C2 #2,3 + 32 Cs #3,0 + 128 Cs #3,1 + 64 C2 #3,2 + 80 C2 #4,o 
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+ 48 C 2 #4,1 + 92 # 6 , + 114 C 4 #o,o,o + 24 Cs #2,0,0 + 48 C2 #2,2,0 + 48 # 2 , 4 ,o 

+ 64 C2 #3,0,0 + 64 C2 #3,1,0 + 64 #3,3,0 + 80 #4,2,0 + 80 #5,0,0 + 80 #5,1,0 

+ 24 ("3 #0,0,0,0 + 48 C2 #2,0,0,0 + 48 #2,2,2,0 + 24 #2,3,0,0 + 48 #2,3,1,0 + 64 #3,1,2,0 

+ 32 #3,2,0,0 + 64 #3,2,1,0 + 64 #4,0,0,0 + 24 #4,1,0,0 + 48 #4,1,1,0 + 92 C2 #0,0,0,0,0 

+ 24 #2,2,0,0,0 + 48 #3,0,0,0,0 + 32# 3) i )0 ,o,o + 36 #2,0,0,0,0,0 + 92 #0,0,0,0,0,0,0 • 

(3.37) 

Finally, at five loops we obtain 

4 5) (x) = - y C 3 3 " 144 a Cs - 204 Cs Ce - 72 C2 Ct - ^ Co - 64 C2 C 3 2 #0 - 80 Cs Cs #0 

- ^ Cs #0 - 384 C3 U #2 - 192 C2 C5 #2 - 288 ( 7 H 2 - 96 Cf #3 - 612 Ce #3 

- 576 C2 Cs #4 - 528 C 5 #4 - 1776 C 4 #5 - 1216 Cs #6 - 1568 C2 #7 - 456 Cs C 4 #0,0 

- 144 C2 C5 #0,0 - 84 C 7 #0,0 - 64 C3 2 #2,0 - 408 Ce #2,0 - 384 C2 Cs #2,2 - 576 C 5 #2,2 

- 768 C 4 #2,3 - 704 Cs #2,4 - 640 C2 #2,5 - 384 C2 Cs #3,0 - 240 C 5 #3,0 - 576 C2 Cs #3,1 

- 864 Cs #3,1 - H52 C 4 #3,2 - 1056 Cs #3,3 - 960 C2 #3,4 - 1656 C 4 #4,o - 1152 C 4 #4,1 

- 1440 Cs #4,2 - H52 C2 #4,3 - 704 Cs #5,0 - 1856 Cs #5,1 - 1216 C2 #5,2 - 1808 C2 #e,o 

- 960 C2 #e,i - 2144 # 8i0 - 48 C 3 2 #0,0,0 - 948 Ce #0,0,0 - 256 C2 Cs #2,0,0 - 160 C 5 #2,0,0 

- 576 C 4 #2,2,0 - 768 Cs #2,2,2 - 384 C2 #2,2,3 - 256 Cs #2,3,0 - 1024 Cs #2,3,1 

- 512 C2 #2,3,2 - 640 C2 #2,4,0 - 384 C2 #2,4,1 - 736 # 2 , 6 ,o - 1368 C 4 #3,0,0 

- 864 C 4 #3,1,0 - H52 Cs #3,1,2 - 576 C2 #3,1,3 - 384 Cs #3,2,0 - 1536 Cs #3,2,1 

- 768 C2 #3,2,2 - 960 C2 #3,3,0 - 576 C2 #3,3,1 - 1104 #3,5,0 - 448 Cs #4,o,o 

- 384 Cs #4,1,0 - 1536 Cs #4,1,1 - 768 C2 #4,1,2 - 1152 C2 #4,2,0 - 576 C2 #4,2,1 

- 1392 #4,4,0 - 1648 C2 #5,0,0 - 1216 C2 #5,1,0 - 384 C2 #5,1,1 - 1648 #5,3,0 

- 1808 #6,2,0 - 2352 # 7 ,o,o - 1568 # 7 ,i, - 368 C2 Cs #0,0,0,0 - 152 Cs #0,0,0,0 

- 912 C4 #2,0,0,0 — 192 C3 #2,2,0,0 — 384 C2 #2,2,2,0 — 384 #2,2,4,0 — 512 C2 #2,3,0,0 

- 512 C2 #2,3,1,0 - 512 #2,3,3,0 - 640 #2,4,2,0 - 640 #2,5,0,0 - 640 #2,5,1,0 

- 288 C3 #3,0,0,0 — 288 C3 #3,1,0,0 — 576 C2 #3,1,2,0 — 576 #3,1,4,0 — 768 C2 #3,2,0,0 

- 768 C2 #3,2,1,0 — 768 #3,2,3,0 — 960 #3,3,2,0 — 960 #3,4,0,0 — 960 #3,4,1,0 

- 1392 C2 #4,0,0,0 — 768 C2 #4,1,0,0 — 768 C2 #4,1,1,0 — 768 #4,1,3,0 — 1152 #4,2,2,0 

- 1184 #4,3,0,0 - 1152 #4,3,1,0 - 1216 #5,1,2,0 - 1376 #5,2,0,0 - 1216 #5,2,1,0 

- 2080 #6,0,0,0 — 1440 #6,1,0,0 — 960 #6,1,1,0 — 2172 C4 #0,0,0,0,0 — 192 C3 #2,0,0,0,0 

- 384 C2 #2,2,0,0,0 — 384 #2,2,2,2,0 — 192 #2,2,3,0,0 — 384 #2,2,3,1,0 — 512 #2,3,1,2,0 

- 256 #2,3,2,0,0 — 512 #2,3,2,1,0 — 512 #2,4,0,0,0 — 192 #2,4,1,0,0 — 384 #2,4,1,1,0 

- 1104 C2 #3,0,0,0,0 — 576 C2 #3,1,0,0,0 — 576 #3,1,2,2,0 — 288 #3,1,3,0,0 — 576 #3,1,3,1,0 

- 768 #3,2,1,2,0 — 384 #3,2,2,0,0 — 768 #3,2,2,1,0 — 768 #3,3,0,0,0 — 288 #3,3,1,0,0 

- 576 #3,3,1,1,0 — 768 #4,1,1,2,0 — 384 #4,1,2,0,0 — 768 #4,1,2,1,0 — 960 #4,2,0,0,0 

- 288 #4,2,1,0,0 — 576 #4,2,1,1,0 — 1728 #5,0,0,0,0 - H36 #5,1,0,0,0 - 192 #5,1,1,0,0 
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— 384 7/5,1,1,1,0 — 368^3 -f^o,o,o,o,o,o _ 736 C2 7/2,o,o,o,o,o — 192 7/2,2,2,0,0,0 — 384 7/2,3,0,0,0,0 

— 2567/2,3,1,0,0,0 — 2887/3,1,2,0,0,0 — 576 7/3,2,0,0,0,0 — 384 7/3,2,1,0,0,0 — 1408 7/4,0,0,0,0,0 

— 576 7/4,1,0,0,0,0 — 384 7/4,1,1,0,0,0 — 2144^2 7/ ,o,o,o,o,o,o — 2887/2,2,0,0,0,0,0 

— 11047/3,0,0,0,0,0,0 — 432 7/3,i,o,o,o,o,o — 736 7/2,0,0,0,0,0,0,0 — 2680 7/ ,o,o,o,o,o,o,o,o • 

(3.38) 

Here all HPLs are understood to have argument x 2 . Note that all indices are positive, in 
other words only the basic indices and 1 appear. This is remarkable, and such a rewriting 
is in general not possible for individual terms contributing to ( |3.34j ) - ( |3.38 ). 

It is very remarkable that within each of the equations ( p.34| ) - fl3.38| ) all terms have 
the same sign, and the common sign is alternating as the loop order increases. In fact, there 
is a sign constraint from the fact that the loop integrals leading to Slo should be positive, 
at least in the Euclidean region < x < 1. Noting that the ladder diagrams appear with a 
factor of (— 1) L per loop order, This implies that (— i) L £l^ is positive for any < x < 1. 
However, the fact that all signs within each of the above expressions are identical seems to 
be a less trivial statement. 

One more check that can be performed on the flff is the limit (f> — > 0, corresponding 
to x — > 1. In this limit the contribution of the ladders to the cusp anomalous dimension 
was derived to all loop orders and to second order in cfr in [10] and reads 

& / 1 -4- ,1k. + 1 \ 



-\lad 



1 - 



2 
1 
4 



24 
7 
512 



16 

+ K 

4387<ft : 
" 497664 



1 + K + 2y/K + 1 

1 50 2 
16 + 288 

21 

+ 0(K 6 ,(t> A 



1 

32 



21 



3456 



5 
256 



2U(f> 2 
20736 



(3.39) 



with k = X/tt 2 . In order to verify this expansion we note that the ladder contribution to 
the cusp anomalous dimension is given by 



ilad 



L>1 



A 



2 ( 8vr 2 



sip , 



and that in the limit we are interested in 

A 



8vr 2 



K . 



(3.40) 



(3.41) 



2(x 2 - 1) 

Taking into account that x = e 1 ^ we expand ( |3.40| ) to second order in </> and find perfect 
agreement with ( |3.39 ) through to five loops. In the next section, we will discuss the limit 
x — > 0. 

3.6 Simplifications in the x — > limit 

The limit x — > is interesting because it connects the velocity-dependent cusp anomalous 
dimension discussed here with the light-like cusp anomalous dimension. 3 

3 Since we have taken the scaling limit we only have a subset of the usual diagrams. However, it is still 
interesting to discuss their behaviour. 
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At four loops, taking the x — > limit of eq. ( 3.37 ) leads to 
(4) {x) x^o 736 lQg7 x + 184** lQg5 x + i6C3 bg4 x + 76^ lQg3 ^ + /32 ^ + ^ lQg2 x 



+ ( "315" + 16 C 3 J log » + Cs + 4vr^5 + 36C 7 ) + 0(x) . (3.42) 

At five loops, we find 

(5) x_,o 2144 9 17152 7 1472 6 2896 5 

fi o ( x ) = ~ TT^T lo § x ~ ~^7F~ & log x - — — C 3 log x — C4 log x 

567 315 45 5 

- (™ C2 Ca + ^ Cs) log 4 x - (64 Cl + 1264 Ce) log 3 x 

- (912 Cs C4 + 288 C 2 Cs + 168 (7) log 2 x 

- (l28 C2 Cl + 160 Cs Cs + ^ Cs) log x 

32 340 
-yCl- 144 C4 C5 - 204 C 3 Ce - 72 C2 C 7 - — Cq + 0(x) . (3.43) 

Finally, at six loops, one obtains 

(6) ^0 339008 n 339008 ^ 9 4288 8 

S2n ' m = log x H C2 log x H C3 log z 

v ; 51975 s 2835 s s 63 s & 

12800 . , , / 34304 , _ 10688 



+ C4 log' x + (^^— C 2 Cs + —jjr- Cs ) log°x 

+ ( ^ Cl + Ce) log 5 x + (5792 Cs C 4 + 1376 C2 C 5 + 528 ( 7 ) log 4 x 

/2944 A o 2432 , 80048 \ , , 
+ ( — C2 Cs + — Cs Cs + — ^— Cs J log 3 x 

+ (128 C| + 3792 C 4 Cs + 7584 Cs Ce + 1152 C2 Ct + 664 &) log 2 x 

8292 

+ (1824 C| C4 + H52 C2 Cs Cs + 336 C| + 672 Cs C7 + C10) log x 



256 , n 2480 

— C2 Cl + 160 Cl Cs + 612 Cs Ce + 432 C 4 C 7 + — 

680 



+ — C 2 C9 + 372Cn+0(x). (3.44) 

It is worth noting that in (|3.42|) - Q3.44 ) certain transcendental constants which cor- 
respond to Multiple Zeta Values [27] having negative indices - such as log(2) or Li4(4) - 
do not appear. This becomes obvious from eqs. (|3.37| ) and ( |3.38|) at four and five loops, 
respectively, but also holds true at six loops. Moreover, eqs. (|3.42|) - ( [3.44 ) contain only 
single zeta values and products thereof. No Multiple Zeta Values of depth 2 or higher 
appear up to six loops, although constants like Cs,3 would be allowed in principle. 

We would like to mention that there is a shortcut for obtaining the asymptotic limit, 
without having to use the algorithm presented above. It suffices to notice that to logarith- 
mic accuracy as x — > 0, we can make the following replacement of the integration kernel 
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appearing e.g. in eq. (|3.6|) , 

1 1 

w' + 1/x w' + X 

Next, rescaling all integration variables by x, we see that one can write the result in the 
small x limit at any loop order in terms of HPLs with indices 0,-1, and argument 1/x. 
The latter can be rewritten in terms of HPLs of argument x, and their small x asymptotic 
behaviour can be made manifest using algorithms implemented in [28]. 



1 



w' + x 



(3.45) 



4. NLO terms in large £ limit 
4.1 Triangle- ladder diagrams (b) 

We now wish to study the sum of the triangle-ladder diagrams shown in Fig. Kb) m a 
similar way to LO. Let F now denote the sum of the diagrams of Figs. |l](a,b), starting with 
1 (as at LO). Then F satisfies the Bethe-Salpeter equation of Fig. [2], with the last term 
omitted. (The last term will be discussed in the following section.) 
Proceeding as at LO, we obtain the differential equation 

8 a d T F(a, t) = Q(a, r) + F(a, r)P(a, r) . (4.1) 

Here the essential new feature is the appearance of Q(a,r). It arises from the first term 
in the second line of the r.h.s. of the equation shown in Fig. |2|. It is given by the one-loop 
integral 

Q(a,T)=c\Xe(^) f^l 1 = c AA^$wf4'4V (4-2) 

J ITT 2 X\ [XX - Z X y{xx ~ Z2) A Zf 2 \zf 2 Z{ 2 ) 

where z\ = e a p >1 and z 2 = —e T q^ are points along the Wilson line, and c = 2/(87r 2 ) 2 . The 
function <&W is known analytically, and we will give a useful form for it later in this section. 
Plugging in the expressions for Zi,z 2 , we have 



cosh(r — a) + cos 4> \ cosh(r — a) + cos 4> ' cosh(r — a) + cos < 
Making the same ansatz as at LO, F = ^2 n e~ nn ^ y2 ^ n (yi,(j)), we obtain 



(4.3) 



.32 A 1 | O 2 (0) 

yi 8n 2 (cosh ?/i + cos 4>) 4 



$ (D ( ., ) . (4.4) 



(cosh yx + cos 4>) \ cosh yi + cos <p ' cosh yi + cos < 

We see that the essential new feature w.r.t. to the LO case is the appearance of an 
inhomogeneous term. It is important to realize that we would like to solve this equation 
to all orders in A, but only to linear order in a = A/A, corresponding to the NLO case. 
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For simplicity of notation, let us abbreviate the potential by — XV and the inhomoge- 
neous term by aXQ. Then we have 



)2 



*( yi A) = aX 2 Q. (4.5) 



Proceeding as in the homogeneous case and setting ^ = e ^V 2 ^ 1 ^ we have 

-d 2 yi r] - -fld yi ri - XVt] = ae +n/2yi X 2 Q . (4.6) 

Recall that at a = 0, this is just the equation for the ladder diagrams, which we already 
solved. We need the solution to order a. We can expand 

V = bladders + OL Va , H = ^ladders + OL , (4.7) 

to obtain, at order a, 

-d$ lVa - ~n ladders d yi r, a - XV 'rj a = e y^^ 2 X 2 Q + ina^dde™ • (4-8) 
As before, Vt is obtained by requiring that ^f'{yi) vanishes at y\ = Therefore we have 

= 2^^7/1^=0. (4.9) 

At order a, this gives 

n a =2a ul (-^—)\ vl=0 . (4.io) 

V bladders/ 

In summary, we have arrived at a differential equation, eq. ( [4.8[ ), together with ( [1.10 ), for 
the contribution of the triangle-ladder diagrams shown in Fig. |l](b). 

We will now explain how to solve these equations to any order in A. First of all, it is 
clear that we can integrate order by oder in A just as we did at LO. The main question is 
whether we can express the resulting wavefunction at each order in terms of the same set 
of iterated integrals as in the previous section. We will now show that this is indeed the 
case, and in fact is true also for a more general class of diagrams. 

The new feature of eq. ( |4.8[ ) is the appearance of Q, so we need to analyze whether 
integrals over Q will be of the same form as at LO. An example will suffice to see that 
this in indeed the case. Consider expanding to order A 2 . Then rf 1 ' 1 {w,x) is given by an 
integral of the form 

— iv W ., u e -V 2 -V (4.11) 



(cosh y + cos </>) \ cosh y\ + cos (f> ' cosh y\ + cos i 

We will now make use of the fact that <&(^ is a function with very special properties. In 
fact, this allows us to immediately make a generalization where 3>W is replaced by <3?( n ). 
This function is given by a beautiful formula [29], 

*W (*,!/) = -==L====& n \x,y), (4.12) 
V(l -x-y) z - 4xy 
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where 



& L) (x,y) = j^pf L L _ // l °S f (z lZ2 ) [U 2L _ f ( Zl ) - U 2L _ f (z 2 )] , (4.13) 

and 

x = ziz 2 , y = (1 - zi)(l - z 2 ) . (4.14) 



Changing variables to w' = e y and x = e 1 ^, eq. ( 4.11 ) becomes, up to a trivial normalization 
factor, 

i ^ ^ (415) 

o w' \w' 2 + 2w'cos4> + 1' w' 2 + 2w'cos(p + I J K 1 



Inspection shows that the variables defined in ( [4.14 ) are given by 



zi = ~ j , z 2 = — l — . (4.16) 
X + w 1 + xw' 

Furthermore, the functions above can be defined using only iterative integrals corresponding 
to symbols z\, z 2 , 1—Zi, l—z 2 . It is easy to verify that the latter factorize over x, w, w+x, 1+ 
wx, and hence are contained in the function class discussed in the previous section. This 
implies that we can again perform all iterated integrals within the set of polylogarithms 
defined by the same integration kernels/symbols as in the homogeneous case, and therefore 
allowing for an algorithmic solution of this problem. 

We note that there is an obvious generalization to a class of diagrams where $W is 
replaced by &( n \ see Appendix A of ref. [10]. The perturbative solution for that class of 
diagrams can be done in the same way as explained above. 

4.2 H-exchange diagrams (c) 

The diagrams with H-exchange of Fig. |l](c) were analyzed in ref. [23] for the quark-antiquark 
potential. It was found that the Bethe-Salpeter equation in that case contains a new term 
of the form 

OO /"OO q 

du dve~^ {u+v) f(u,v; yi )^( yi -u + v), (4.17) 
o Jo 

so that one has a linear integro-differential equation for \F Their analysis can be adapted 
to the present case of general (j), with / now depending on <f>. 

Although such an equation may seem complicated, it simplifies considerably solving 
it in the small a = A/A limit. The reason is that the kernel, the .ff-exchange diagram is 
already of order a, so that we only need the wavefunction at order a . In other words, 
the problem reduces to a differential equation for the wavefunction at order a, with an 
inhomogeneous term. This is exactly the case we studied in the previous section. 

Having said this, the main difficulty lies in the computation of the H insertion, and in 
integrating it when iteratively solving for the wavefunction. From the discussion above it 
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is clear that we need to understand how to carry out the H-shaped and similar integrals. 
Let us therefore start with the basic three-loop integral, which has one H-exchange, and 
no additional rungs. It is given by 

poo rS2 roo i-t2 

/ ds 2 / d8! / dt 2 / dt 1 f(-s 1 ^,-s 2P ^,t 1 q^t 2 q^. (4.18) 
Jo Jo Jo Jo 

Note that strictly speaking we should introduce IR and UV regulators for this integral, but 
since we are only interested in extracting the overall divergence, the details of the cutoffs 
are not very important. For the same reason, the H-shaped subintegral can be defined in 
exactly four dimensions, 



f(x 1 ,x 2 ,x 3 ,x 4: ) =(<9i + <9 4 ) 2 h(xi, x 2 ; x 3 , x±) , (4.19) 
d 4 x^d 4 XQ 1 

(/i^77-2'\2 />! 2 ,-y> 2 /y«2 ^y>2 
in ) x l5- t '25 x 36 x 46 x 56 



h(x u x 2 ; x 3 , x 4 ) = / , 2 . 2 2 2 2 2 2 • (4.20) 



Eq. ( 4.19 ) defines the function /. Although this is a two- loop integral, / reduces to one-loop 
integrals thanks to differential equations it satisfies. We review these differential equations 
in Appendix A. Remarkably, they allow us to express / in terms of the one-loop function 
the same function that appeared already in integral class (b). Explicitly, we have 

/ = xUx\ 2 + xj 3 - xl) (4. 40 + 4.(4. + Aa - 4d ( 4) 

\X 13 X13/ \ x 24 x 24/ 

_L T 2 C„2 , 2 _ 2 x *(1) ( 4± £o4 \ , 2 / 2 , 2 _ 2 x ^(1) / £34 x 23 A 

V x 13 ^13 / \ x 24 x 24/ 

, / 2 2 2 2 _ 2 2 \ ^(1) ( X 12 X 34 x 14 X 23 \ (a 91 \ 

+ l x 13 x 24 x 14 x 23 X 12 X 34J ^ I 2 2 ' 2 2 ) ' l^- zi / 

\ X 13 X 24 X 13 X 24/ 

where / = (^i2 x i3 x 24 x 34)/- This formula will be very convenient when discussing the 
strong coupling limit. 

After this digression on h, we can proceed to extract the overall divergence and compute 
the H-exchange integral. Changing variables according to s\ = x\s 2 ,t\ = x 2 t 2 , and s 2 = 
zp,t 2 = pz, where z = 1 — z, and using that h scales as l/x 4 , we find 

[°°^H^, (4.22) 
Jo P 

where 

= [ dzdx 1 dx 2 f{-x 1 zp tI ,-zp fl ;x 2 zq fl ,zq^). (4.23) 
Jo 

Note that by assumption is finite (i.e. the original integral only had an overall UV 
divergence). However, when carrying out the integration in ( [4.23| ), care is required, be- 
cause the finiteness is not necessarily true for individual terms appearing in ( 4.21| ). This 
small problem can be avoided by introducing an auxiliary regulator. With the above 
parametrization, we have 

= -xizpf 1 , x% = -zpf 1 , = x 2 zq» , x% = zq^ 1 , (4.24) 



- 18 - 



and using p 2 = q 2 = 1 ,p ■ q = cos <fi, we have 



x\ 2 =x\z 2 , X23 = Z 2 + x\z 2 + X2ZZ1 cos <f> , (4-25) 

and so on. 

In summary, we found a finite parameter integral, where the number of integrations 
equals the expected degree of the function. Just as for integral class (b), higher orders 
can be obtained by iteration. However, it is not yet clear that the same class of functions 
will be sufficient to evaluate these integrals. Explicit results at three loops motivate that 
it might be. We leave this question for future work, and close this section by remarking 



that formula (4.13) will certainly be very useful when trying to evaluate this integral and 



similar integrals appearing in the iterative solution. 
5. Strong coupling limit at LO and NLO 

Here we discuss the strong coupling limit of the Bether-Salpeter equations. In this limit, 
the calculation of the ground state energy becomes almost trivial. It is straightforward to 
extend the analysis of ref. [23], which was done in the anti-parallel lines limit <j) — > it, to 
any angle. 

5.1 Strong coupling limit of Bethe-Salpeter equation 

Let us start by discussing diagrams of type (c). First of all, we notice that as in [23], 
the Bethe-Salpeter equation for this class of diagrams simplifies dramatically in the strong 
coupling limit. The reason is that for fio ~ Va S> 1, the region of small u,v will give the 
dominant contribution to the integral in eq. ( [4.17 ). This implies that the wavefunction 



fy(yi) can be pulled out of the integral, with the coefficient being an effective potential. 
This argument also works for the angle-dependent cusp Wilson loop. 

We therefore need to compute the effective potential for general angles. Although the 
function h is not known analytically, its derivative / is known, as we saw in the previous 
section. 

We need the function f(xi,X2] X3, X4) in the limit where x\ — > X4 and X2 — > X3. Let us 
parametrize this limit by x 2 4 = u 2 ,x 2 3 = v 2 ,x 2 2 = x 24 = x i3 = x 34 = 2coshyi + 2cos(/>, 
with u,v small. Plugging these values into eq. ([4.21 ), it turns out we only need the 
following limit of 



$«(!,£) = _i oge + 2 + 0(e) . (5.1) 



Using this limit, we obtain 



, , . u 2 log u + v 2 log v 2 2 

t — > —4- ~T"7r + 0(u , v ) . (5.2) 

J (2 cosh y x + 2 cos c/>) 3 v ' ; y ' 

We see that this is a generalization of eq. (5.2) of [23] to general angles. One could also 
use eq. ( [4.21| ) to compute higher order terms in the expansion. 
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This means that the correct effective potential for the general angle case is obtained by 
replacing each [x 2 + 1) terms in (5.3) of [23] by (2coshyi + 2cos</>) for the cusped Wilson 
loop. Then we have a Schrodinger equation 

n 2 ' 



-d 2 yi + v eS ( yi )+ 4 



*(yi) = o 

where the correction to the effective potential comes from the integral 



V, 



efflAA 2 



dudv e ~^ u+v) f(u,v) . 



Explicitly, we have 
V eS = 



A 



+ 



AA 2 logO 



(5.3) 



(5.4) 



(5.5) 



47r 2 (2coshyi + 2cos0) 2vr 6 C7 4 (2 cosh y x + 2cos</>) 3 ' 

At strong coupling, we can focus on A S> <C 1, with A(yi) 1//4 fixed. In that regime the 
leading term of the Schrodinger equation is 



v cS ( yi = 0) + 

From this we obtain for the ground state energy. 



O 2 

5 i 



0. 



r (a)+(c) 



-On 



2"7T COS ■ 



1A, A m 
2l l0g A + ° 



(5.6) 



(5.7) 



Here the superscript indicates that this is the contribution from the integrals shown in 
Figs. 0(a), (c). 

Let us now discuss the integrals of Fig. |l](b). Here we obtained a Schrodinger equation 
with an inhomogeneous term (note that there a = A/ A) that is not multiplied by the wave 
function. The latter fact suggests to us that the contribution of this class of diagrams at 
strong coupling will not be given by an exponential factor of the type seen for integral 
classes (a) and (c). If one assumes the absence of contributions of integral class (b) at 
strong coupling, as we will do in the following, then ( |5.7D is the full answer at LO and NLO 
in the scaling limit. 

Let us now compare this against the corresponding quantity computed in string theory. 
5.2 Scaling limit of the string theory result 

The leading term (and first subleading term as well) in the \j\f\ expansion at strong 
coupling has been computed using string theory in ref. [7]. It is straightforward to expand 
their result in the large A limit that we are interested in. For the LO, this was already 
done in ref. [10]. 

It is easy to take the scaling limit of the formula for T given in ref. [7]. The details of 
this calculation are presented in Appendix B. We find 



2tt cos I 



1 A A , n 
2l bg A + ° 



(5.8) 
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As a consistency check, we can take the limit <j> = ir — 8, S — >■ 0, where we expect to find the 
quark-antiquark potential V. More precisely, T ~ —1/SV, and indeed we find agreement 
with eq. (5.4) of [23]. 

Let us compare eq. Q5.8Q to the diagram calculation performed in the previous sub- 
section. Comparing to eq. ( |5.7|) , we find perfect agreement. Recall that in principle there 
could also be a contribution from integrals of type (b) not accounted for in eq. ( |5.7| ) , but 
we argued that this is not the case based on the structure of the Bethe-Salpeter equation 
for these integrals. Under this assumption, we see that there is a perfect match between the 
field theory calculation in the scaling limit, and the string theory calculation. As pointed 
out in [10], this agreement was not guaranteed due to potential order of limits issues. 

6. Discussion and conclusion 

In this paper we further studied the scaling limit of the cusp anomalous dimension intro- 
duced in [10], in several ways. 

In the first part of the paper, working at LO we showed that the perturbative solution 
at weak coupling can be expressed at any loop order in terms of harmonic polylogarithms, 
and outlined a corresponding algorithm. As illustration, we reproduced the three-loop 
result of [10] and computed the four- and five-loop results, which are new. We also provide 
a shortcut for obtaining the x — > asymptotics, which corresponds to the light-like limit 
of the edges of the Wilson loop. 

We observed interesting features of these results. We find that, at least up to five 
loops, they can be written in terms of a reduced class of harmonic polylogarithms, with 
indices and 1 only, when choosing x 2 as argument (this feature was already noted in [10] 
up to three loops.). Moreover, in the x — > limit, up to six loops, we find that the resulting 
asymptotic expansion can be expressed in terms of linear combinations of products of single 
zeta values only. Other constants such as log(2), or multiple zeta values of higher depth 
were not needed. This is especially interesting in the context of the BES equation for the 
closely related light-like cusp anomalous dimension, which has the same property [30]. 

It would be very interesting if one could prove these properties. Such a proof would 
likely shed more light on the structure of the cusp anomalous dimension. 

In the second part of the paper, we extended the analysis of [10] to NLO order. The 
new feature of the equations is the appearance of an inhomogeneous term. (A similar 
analysis was recently done for the quark-antiquark potential [23]). This term does not 
alter the perturbative solution, however, and we were able to apply the same strategy as at 
LO. We showed how to compute these contributions systematically in perturbation theory. 
For one class of integrals, we provided an algorithmic solution at any loop order in terms of 
harmonic polylogarithms. For the second class of integrals, we showed how to obtain the 
solution in terms of iterated integrals of simple functions. We left the question of whether 
the latter can be expressed in terms of HPLs for future work. 

Finally, we discussed the strong coupling limit of the equations. We computed the 
logarithmically enhanced terms at NLO, and found agreement between the field theory 
and the string theory calculation. This generalizes the calculation of [23] to any angle. 
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Using our formulas, it should be possible to compute the non-logarithmic terms at NLO 
as well. We leave this for future work as well. 

It would be interesting to extend the analysis of [10], in the zero angle case to NLO. 
In that case, the Schrodinger potential becomes the exactly solvable Poschl- Teller poten- 
tial. At NLO order, one has the same homogeneous part of the equation, and a new 
inhomogeneous term, as discussed in the present paper. 

Our approach also suggests a general strategy for the computation of the cusp anoma- 
lous dimension, or related quantities. At a given loop order, there are two sets of contri- 
butions. First, there are a number of integrals that have an overall UV divergence. These 
diagrams are the "seed" of the Bethe-Salpeter equations and have to be computed. They 
correspond to the most complicated part of the calculation. However, the fact that they 
have no subdivergences allows one to extract the overall divergence easily, so that one is left 
with the calculation of a finite quantity. The latter is sometimes related to four-dimensional 
integrals. This observation allowed for example for the computation of an infinite class of 
integrals contributing to the cusp anomalous dimension in ref. [10]. Second, there are di- 
agrams that do have subdivergences. For these contributions, the resummation technique 
via the Bethe-Salpeter equation is very useful, as it automatically takes into account the 
non-Abelian exponentiation. Although these contributions typically give the most compli- 
cated contributions as far as the functions involved are concerned [10], the latter have their 
origin in simple iterations of diagrams of the first type. 

Although our analysis did not rely on the planar limit, non-planar contributions to 
r cU sp appear only at four loops, or at higher subleading terms in the scaling limit. It would 
be very interesting to compute the first non-planar corrections. We expect that many 
observations about the calculation of loop integrals, especially the comments for extracting 
overall divergences and using four-dimensional integrals, will be useful in related problems 
as well, e.g. as the non-planar integrals discussed in ref. [31]. 

Our approach can also be extended beyond the NLO. We remark that this does not 
require any Feynman graph calculations, as the integrand for the planar Wilson loop can be 
obtained from a soft limit of the integrand of a four-particle scattering amplitude [10,11]. 
The latter can be obtained through on-shell recursion relations in principle to any loop 
order. We give examples of this procedure in Appendix C. We hope that this all-loop 
knowledge of the Wilson loop integrand gives a good starting point for analyzing further 
the properties that we have observed in this paper. 

Finally, the scaling limit discussed here might be useful for simplifying the TBA equa- 
tions of refs. [16,17]. It would also be interesting if those equation could shed light on some 
of the observations about the perturbative properties of r cusp that we have made here. 
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A. Differential equations for two-loop integral h 



Here we give the differential equations for the two-loop integral h of eq. ( (4.20| ), 

The non-trivial differential equation has been written down e.g. in eq. (A. 7) of ref. [32]. 
The differential operator we have can be related to the one of eq. (A. 7) of [32] by using 
translational invariance Q^i=i ®i = 0), up to trivial pieces proportional to where we 
can use the Laplace equation. Note that all occurring terms can be written in terms of 
thanks to eq. (A.5) below. 
More explicitly, we have 



-2(di + 3 4 ) 2 = -(di - d 2 ) ■ {d 2 - d 3 ) - (8 1 + d 2 f + Di + D 2 + D 3 + D 4 
We have, using eq. (A. 7) of ref. [32], up to overall factors, 



(A.l) 



(d l - d 2 ) ■ (d 2 - a 3 )h 
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x 12 x 34 



(X 13 X 24 - X 14 X 23 )Xi234 + (X U - X 13 )Xi 34 

- (x 24 - ^23)^234 + ( x l 2 - a?li)X 3 i2 - (x 42 - x 41 )X 4 i2 . (A. 2) 
and, from the Laplace equation, 

[-{di + d 2 f + Ql + D 2 + D 3 + D 4 ] h = - X l23A + -^(X 1U + X 2U ) + -^(X 123 + X 124 ) , 
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_$M ±12^^23 
x 13 \ x 13 x 13 



(A.3) 

(A.4) 
(A.5) 



Combining differential equations ( [A. 2; ) and ( |A.3| ), and plugging in ( |A.4 ) and ( |A.5 ), we find 



eq. (4.21) given in the main text. 



B. Details of the strong coupling calculation 

Here we show the details of the expansion of the string theory answer for V in the scaling 
limit. The result of [7] for V is parametrized by two parameters p and q = —ir , r > 0, 
which are implicitly defined through the angles <f> and 9, in the following way 

2bq 



y/W+p 2 



K(k 2 

2 



-IT 
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b^+p 2 



n 



b A +p- 



K(k 2 



where 



(,2^2 



q 2 ) 2 + 4p< 



b 2 (b 2 -p 2 ) 
b A +p 2 



(B.l) 
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(B.3) 
(B.4) 
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In terms of these variables, we have 



_V\ 2^/b r +p 2 Ub 2 + l)p 2 
2vr bp [ b A +p 2 

where E, K and II are complete elliptic integrals, 

1 



K{k 2 ) - E{k 2 ) 



E{k 2 ) 
K(k 2 ) 
n(a 2 ,6 2 ) 



dt- 



\fl — k 2 sin 2 1 
dt\J\ - k 2 sin 2 t, 



dt 



(1 - a 2 sin 2 t) \A - b 2 sin 2 1 



(B.5) 

(B.6) 
(B.7) 
(B.8) 



The scaling limit i9 3> 1 is reached by letting p — > 0. We see that we require the leading 
and subleading divergences of II in the limit where 
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7T — 2 arcsin 



lim 11(1 — ae, 1 — eb) - _ 
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log(e) 
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In this way, we obtain, at leading order in \/A S> 1 
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2l (1+r 2 ) 
1+p logp ^ 4 
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(B.10) 



We now convert r and p to their expressions in terms of 9 and (f). So we need the 
expansions of the latter to the necessary order in p. We hnd 
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P Vl + r 2 



+ plogp 



3 + r 2 



r 2 \/l + r : 



+ 0(p). 



=2 arcsin ■ 



1 



1 



p 2 logp + 0(p 2 ) , 



We can see that this is in agreement with equation (|5.8|) quoted in the main text. 



(B.ll) 
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C. Relation between integrals for four-particle scattering amplitude and 
cusp anomalous dimension 

In Figs. |3| and |4| we illustrate the relation between the integrals contributing to the four- 
particle scattering amplitude (odd lines) and the integrals contributing to the cusp anoma- 
lous dimension (even lines) to three loops. The integrals occurring at one and two loops 
are shown in Fig. ||. The three-loop integrals are shown in Fig. |j. The reason for this 
relation [10,11] is essentially exact dual conformal symmetry [15], together with the fact 
that massive scattering amplitudes in soft limits are related to Wilson loops [33]. 
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Figure 4: Relation between integrals of the four-point amplitude (first and third line) and Wilson 
line integrals (second and fourth line) at three loops. Arrows denote internal numerator factors 
(jy* + q^ 1 ) 2 , where and q^ are the momenta flowing along the lines with arrows. 
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